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In order to increase the accuracy of the linearized augmented plane wave method (LAPW) we 
present a new approach where the plane wave basis function is augmented by two different atomic 
radial components constructed at two different linearization energies corresponding to two different 
electron bands (or energy windows). We demonstrate that this case can be reduced to the standard 
treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions 
of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show 
that the task is closely related with the problem of extended core states which is currently solved 
by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, 
the number of supplemented basis functions in our approach is doubled, which opens up a new 
channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new sup¬ 
plemented basis functions absent in the LAPW+LO treatment is closely related with the existence 
of the ui— component in the canonical LAPW method. We discuss properties of additional tight 
binding basis functions and apply the extended basis set for computation of electron energy bands of 
lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We 
demonstrate that the new treatment gives lower total energies in comparison with both canonical 
LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor 
LAPW basis sets. 


I. INTRODUCTION 

The choice of a basis set which may first appear as 
“the black art” [lj is being constantly debated within the 
quantum chemistry community. It is well known that 
in molecular calculations there are two main groups of 
molecular basis sets introduced by Pople and collabora¬ 
tors (see lj and references therein) and more recently by 
Dunning (2]. Both groups supply us with the whole hier¬ 
archy of basis sets, where at each step we can enrich the 
main set by polarization functions of high orbital type or 
by diffuse functions. For example, one has to add polar¬ 
ization functions if polarization effects are expected to be 
important, or diffuse functions if we want to refine the 
description of extended molecular states. Not surpris¬ 
ingly, the actual choice of a basis set depends on the task 
to be solved and is considered as a difficult problem. For 
heavy and laborous calculations the choice of a basis set 
is crucial since on one hand we want to obtain a reliable 
result and on the other hand minimize the computer time 
to achieve the goal. 

In contrast to this complicated hierarchy of molecular 
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basis sets, the choice of bases in electronic band structure 
calculations and here we imply mainly the linear aug¬ 
mented plane wave (LAPW) method [3H3 , seems rather 
simple. The number of augmented plane waves is com¬ 
monly determined by the parameter RMTK max , where 
Rmt is the smallest muffin-tin (MT) radius, and K max 
is the maximal value of the plane wave vector. The K ma x 
value implies that the kinetic energy cut off is K^ nax /2 
(in atomic units). However, in practice some band struc¬ 
ture calculations can not be carried out without so called 
local orbitals (LO) [g]. Such situations occur in systems 
with semicore electron states which cannot be fully con¬ 
fined within the MT-spheres. The problem of extended 
core states and its relation to our approach is discussed 
in detail in the next section. What concerns us here is 
that the introduction of local orbitals represents an ex¬ 
tension of the canonical LAPW basis set albeit the form 
of new basis states (local orbitals) is very different from 
the standard augmented plane wave basis function. The 
LAPW+LO scheme proposed by Singh Q, j|} is prac¬ 
tical, but the way it has been introduced is not fully 
satisfactory. The form of the local orbital functions is 
not derived from a general approach, and arguments for 
adding LO basis states are purely variational. 

In the present paper we show that the appearance of 
supplemented basis functions can be understood as a re- 
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suit of refinement of the LAPW band scheme when in 
an effort to increase its accuracy we use two linearization 
energies (corresponding to two electron bands). We will 
demonstrate that new basis states are of two types. The 
first group gives the local orbitals in the form suggested 
by Sing Q. however, the basis functions of second type 
have different form which is not used in the LAPW+LO 
method. Therefore, the canonical LAPW basis set and 
also the LAPW+LO basis set can be extended to a more 
complete basis set. New basis functions and consequences 
of their introduction are closely examined in the present 
work. 

The paper is organized as follows. We start with re¬ 
visiting the problem of extended core states which gives 
rise to the LAPW+LO scheme and formulate our ini¬ 
tial statement for the refinement of the LAPW method, 
Sec. nn In Sec. m we present our method which results 
in adding tight binding basis functions to the canonical 
LAPW basis set. In Sec. [TV] we apply the method to 
electron band structure calculations of the face centered 
and body centered lattice of lanthanum and the hexago¬ 
nal close packed lattice of cadmium. Our conclusions are 
summarized in Sec. El 


II. THE PROBLEM OF EXTENDED CORE 
STATES 


The linear augmented plane wave (LAPW) method [j§- 
H] is probably the most precise method for electronic 
band structure calculations and is widely used for the 
calculation of materials properties (B|. 

In the LAPW method space is partitioned in 

the region inside the nonoverlapping muffin-tin ( MT) 
spheres and the interstitial region I. The basis functions 
<l>n{k, R) where n = 1,2,..., N are given by 


(j) n {ki R^) — 


V “ 1/2 exp(i(k + K n )R), 
£ Et)Yi, m {r), 


with radial parts 


Re I 

R e MT{a) 
(1) 


E t ) = A?£ t H (r, E t ) + u t (r, E l ). (2) 

Here the index a refers to the type of atom (or MT- 
sphere) in the unit cell, the radius r is counted from 
the center R a of the sphere a (i.e. r = R — R a ), V is 
the volume of the unit cell. Radial functions iq )m (r, Ei) 
are solutions of the Schrodinger equation in the spher¬ 
ically averaged crystal potential computed at the lin¬ 
earization energy Ei, and ui(r,Ei) is the derivative of 
ui t m with respect to E at Ei. The coefficients and 
.B" m are found from the condition that the basis func¬ 
tion (f> n is continuous with continuous derivative at the 
sphere boundary, r = R^t (Rmt is th e radius of the 
MT -sphere a). In the following for compactness we omit 
the index a and restore it when needed. Linearization 


energies E[ are chosen close to average values of corre¬ 
sponding band energies or to the Fermi level. The ex¬ 
tended electron basis states defined by Eq. © as a rule 
are orthogonal to the core states. This is a consequence 
of the relation 


cRmt 


U v (r) U c (r)r 2 dr = 


d2 

n M T 


U c (R MT ) dUv ^ r MT) - U v (R MT ) 


2(E C - E v ) 
dU c (RMT)\ 


dr ) ’ 


( 3 ) 


applied to a core state with orbital quantum numbers 
l,m and the radial wave function U c [r ), and a partial 
radial function of valence state, U v (jr) = 7 Ei), Eq. 
(ED, with the same angular dependence. Notice that Eq. 
© ensures the orthogonality between the extended and 
core states if two conditions at the sphere boundary are 
satisfied for each of the core states, 


Uc(Rmt) = 0, 
8U c (Rmt) 
dr 


(4a) 

(4b) 


Although these conditions are met for a great number of 
cases, they are violated for so called semicore states that 
are not fully contained in the muffin-tin sphere 0-i- 
Semicore states leaking out of the MT-regions should 
be treated as extended states. This in turn requires that 
the linearization energy Ei is chosen near the energy of 
the semicore level, Ei « E c , because the LAPW basis de¬ 
scribes only states near Ei well. However, as E c is quite 
far from the Fermi energy Ep and the valence band en¬ 
ergy, the choice Ei ss E c inevitably gives poor descrip¬ 
tion for partial l = l c valence states. On the other hand, 
the option Ei = E v is not satisfactory for the semicore 
states situated substantially deeper in energy. As dis¬ 
cussed in Refs. [7 : , jj jlll - llll there is no simple solution 
to this dilemma. Even worse, in many cases the attempt 
to use a single value of Ei for both valence and semi¬ 
core states leads to the appearance of so called “ghost 
bands” Hi giving false band energy positions. As a 
remedy one can divide the energy spectrum in two win¬ 
dows (energy panels) and use two different sets of Ei for 
calculations of semicore and valence states, respectively 
[sj. This technique however is also not fully satisfactory 
because now there is no single Hamiltonian matrix for the 
problem and strict orthogonality between electron states 
belonging to different energy windows is not guaranteed. 
Ideally, in the MT-region there should be two different 
types of radial components with the same angular depen¬ 
dence l = l c . That is, in Eq. © 


K,m (r) = (r) + n\2 n (r), (5a) 

where 

El 1) )=A^ n u[ 1 \r,EP) + B ^ n u , i 1 \r, E l 1) ) 

(5b) 
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refers to the semicore states with Ej 1 ^ = E c , and 


/ jp( 2),n/ -pp(2)\ (2)/ rp(2)\ n(2),n *(2)/ zt>(2)\ 

] \m W W )-\m U l V’h ) + B l,m U l W A, ) 

(5c) 


( 2 ) 

refers to the valence states with E\ = E v ss Ep. Both 
states, Eqs. (IHbl) and (l5cj) . should merge to a single 
l c —wave component of the plane wave 


<f> n (k, R) = V 1 / 2 exp(i(k + K n )R) (5d) 


at the surface of MT sphere. Now, however the boundary 
problem becomes ill-defined, because for the l = l c com¬ 
ponent there are four coefficients, Ap^’™, S^’", A^’” 
and Dpi', ' 1 for only two boundary conditions. 

In Ref. Q Singh has proposed to increase the number 
of boundary conditions to four by matching the value of 
the basis function and its first three radial derivatives at 
the sphere surface. This gives rise to super-linearized 
APW method denoted as SLAPW-4 Q because four 
functions, four coefficients and four boundary conditions 
are involved. In a simpler super-linearized modification 
called SLAPW-3 Q the first radial part 711 [p”, Eq. (IFbl) . 

/Q\ 

is supplemented by only one function u) ' (instead of 
Eq. ([5c])). The three coefficients (A^’", B^ n , 

App") are determined by requiring continuity of the ba¬ 
sis function and its two derivatives. 

In comparison with standard LAPW method in both 
SLAPW modifications there are additional requirements 
for the plane wave convergence. Indeed, the plane wave 
expansion of the interstitial region must converge either 
to the correct second and third derivative (SLAPW-4) 
or to the second derivative (SLAPW-3). Because of that 
more plane waves are needed to satisfy these additional 
requirements, the plane wave energy cutoff parameter 
should be increased and calculations become much more 
costly 0,11- 

To circumvent the problem and improve the LAPW 
efficiency Singh put forward a third approach based on 
local orbitals (LAPW+LO) 0. In the LAPW+LO ap¬ 
proach the same three radial functions as in SLAPW-3 
are used (i.e. u^\ u\ 1 ' > and u | 2 p, but the coefficient of 
U; 2 '* is fixed (say, A^’™ = 1) and the two remaining co¬ 
efficients (A^’”, B^ n ) are found from the conditions 
that the local orbital goes to zero with zero derivative at 
the sphere boundary. Nowadays, LAPW+LO is widely 
used for band structure calculations of solids with semi¬ 
core states Si- However, conceptually the LAPW+LO 
method is understood as a procedure giving additional 
variational freedom through an increase of the number 
of basis functions. It is not clear why additional ba¬ 
sis functions should include these particular components 
(i.e. u^\ u[ 1 ' > and u| 2 ^). The proposed zero boundary 
conditions for local functions are not derived from a gen¬ 
eral physical statement. 

Inspired by the LAPW+LO method 0 in the present 
study we formulate a more general approach to the prob¬ 


lem. Unlike the LAPW+LO approach which uses varia¬ 
tional arguments for its foundation, we will derive sup¬ 
plemented basis states from the initial requirement that 
two different radial functions (i.e. 72.^’ 71 and T^pp", Eqs. 
(l5bl) . (l5cl)l. having the same angular part merge in a sin¬ 
gle plane wave function <j> n , Eq. (I5dl) . in the interstitial 
region. Unlike SLAPW-4 or SLAPW-3 we retain only 
two joining conditions across the MT-sphere boundary. 
As a result, we will obtain two types of supplementary 
tight-binding basis functions (see Eqs. (I15bl) and (I15cl) 
below), satisfying Bloch’s theorem. 


III. DESCRIPTION OF THE METHOD 


As discussed in Sec. El in the case of semicore states 
we have two types of radial solutions in the MT-region 
with the same angular dependence Y^ m (r) but differ¬ 
ent linearization energies e\^ and E^: ^[^’"(r, Ef 1 ^), 

Eq. (t5bl) . and ^[^’"(r, e[ 2 ' > ), Eq. (l5cl) . One of the ra¬ 
dial functions can refer to extended states, i.e. R e (r) = 
K^ n (r, E^), while the other can refer to supplemen¬ 
tary angular states R a (r) = 7£p2i’"(r, E[ 2 ^). As we will 
see later in Sec. El in practice we describe the semicore 
states as extended states with Ei = E core and valence 
states with the same l as supplementary states for which 
Ei = E v . (For metals one can take Ei = E v « Ep.) 
Since in the interstitial I -region both types of solutions 
are represented by the plane wave function 4> n (k. R), Eq. 
(I5dl) . they become indistinguishable there. In the LAPW 
method there are two matching conditions (for the func¬ 
tion and its derivative) on the sphere boundary. There¬ 
fore, in our case we have 


A e u e + B e ii e + A s u s + B s ii s — 

An r A -» -» 

yii(UMr)C(W (6a) 

A e u e + B e u e + A s u s + B s u s = 

Anr s. -» -» 

—7=i l ji{k n R M T)Yi* m {k n ) e“. (6b) 


Here we adopt short notations A e = App", A s = App", 
u e = u[ 1) (R M t,e [ 1) ), u s = uf , \R M T,E\ 2) ), u' e = 

du < j 1 \R M T,El 1 ' > )/dr, it' = duf‘\R M T,E < f , ' l )/dr , and 
have used the Rayleigh expansion of the plane wave <f> n 
on the sphere surface. Since there are four coefficients 
(A e , B ei A s and B s ) and only two equations, it is clear 
that the general solution to Eqs. (l6al) . (Ihbl) . forms a two 
dimensional linear space with two linear independent ba¬ 
sis vectors. 

Further, introducing standard LAPW quantities a e = 
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ap ),n , b e = b[ 1 ' >,n , where 

A ■jr A —» —* 

A e = -j=i l R 2 MT Y* m (k n ) e lknRa a e , (7a) 

/j ' 77 * A —* —* 

b e . (7b) 

and analogous relations for a s , 6 S , we rewrite Eqs. (l6al) . 
(Iffbl) as 

a e u e + b e u e A a s u s A b s u s = ji(k n R MT )-^%—, (8a) 

R mt 

a e u' e A beii'e A a s u' s A b s u' s = j[ (k n R M T)-Jr~ ■ (8b) 

ri MT 

Notice that the standard LAPW solution for a e = a° 
and b e = b° e without supplementary states, i.e. when 
a s = 0, b s = 0, can be found from the following system 

a° e u e + b° e u e = ji ( k n R M T)-jJr ~, (9a) 

ri M T 

o° e u' e + b °eK = (9b) 

U MT 

Defining auxiliary quantities t a and t b 

t a =a e - a°, (10a) 

tb = b e - b° e . (10b) 

and subtracting Eq. (l9all from Eq. (l8al) . and Eq. (t9bl) 
from Eq. (I8bl) we arrive at 

t a u e A t b u e A a s u s A b s ii s = 0, (11a) 

t a u e A tbU e A a s u s A b s u s = 0. (lib) 

The solution to Eqs. (HUD, EB can be found from 
the following two systems, 


f ds,iu s T b s ^ii s — u e 
\ a s ,iu' s A b S}1 u' s = -u'e ’ 

f a s , 2 U s A b S ' 2 Us = —Ue 
\ a s - 2 u' s A b s . 2 u' s = —u'e 


(12a) 

(12b) 


Solutions to the systems (ll2aH and (Il2bl) are quoted ex¬ 
plicitly in Appendix lAl Eqs. (IAlal) - (IA2bl) . Having found 
a s ,i and b s ,i ( i = 1,2), we write the general solution to 
Eqs. HUD, (IFEl) as 

a e — d e A ta , (13a) 

b e = b° e + t b , (13b) 

d s — t a Us,l 4“ t b Us, 2 , (13c) 

b s = t a 6 Si i + tb & s ,2, (13d) 


where t a and t b are arbitrary numbers. The full radial 
component TZf m (r) of the basis function inside the MT- 
sphere a , Eq. (I5al) . is written as 

R'l,m(j') ^ U’e Ue A b e lie A ta {ll e A U S; 1 U s A b s \ il s ) 

+t b (u e A U S) 2 u s A &S, 2 u s )- (14) 


(Here notations u e = u e {r), u s = u s (r) etc. refer to 
radial functions.) 

Notice that since the coefficients t a and t b are arbi¬ 
trary, they should be found from the standard variational 
procedure by requiring the minimization of the LAPW 
ground state energy. Furthermore, the form (fl4l) sug¬ 
gests considering three linear independent radial parts 
(i.e. Bf m (r), i?f)^(r), -Rf’^(r)) instead of the single func¬ 
tion + *0 + tb Ex P 11(> 

itly, 


R tm( r ) = Cn e lknRa (a° e u e A b° e lie), (15a) 

= C n e %knRa (u e A a s ,i u s A 6 s ,i u s ), (15b) 
R a i)m( r ) = Cn e l ^ nR °‘ (u e A a s ,2 U s A 6 S , 2 u a )- (15c) 


Here u e = u e {r ), u s = u s (r) etc. are corresponding radial 
functions and 


C n = -j=i l R 2 MT Y l * m (k n ). (16) 

The first function, Eq. (I15al) . is in fact the standard ra¬ 
dial part of the l— type, Rf m = Eq. 0, entering 

the usual LAPW basis function R), Eq. 0. Its 

coefficients a° e and b° e are given by the LAPW boundary 
relations, Eqs. (f9a|) . (|9bD - Two other functions however 
are very different from <j) n (k, R) and should be included 
to the LAPW basis set as extra basis states, 

<t> s>i {k,R) = Yr{r)R s l 'j m {r), (17) 

where i = 1,2. The important thing is that their coef¬ 
ficients a Sj i, 6 Sj 1 are found from Eq. (I12al) . while coeffi¬ 
cients a Sj 2 , b Si 2 from Eq. (Il2bl) . In respect to two new 
functions Tq’A Eqs. (I12all . (H2bl) impose the following 
boundary conditions 


K» = 0 , 

(18a) 

dR^Ur) 

dr 

(18b) 


These relations have a simple interpretation: new sup¬ 
plementary basis functions 4 > s ,i(k 1 R) are required to be 
orthogonal to the standard LAPW radial functions. We 
want to stress here, that the conditions (IISail and (Il8bl) 
are not assumed or introduced at our will. They are de¬ 
rived from the initial equations hid, m [or equivalently 
from Eqs. (|8aP . d8bP ] and are used to obtain the general 
solution, Eqs. (I13all - (ll3dl) . 

From Eq. (I15bl) , (H5cD and (1161) it follows that the sup¬ 
plementary basis states in principle depend on the in¬ 
dex n, i.e. 4> s ,i(k, R) = cK,i{k n , R)- (We recall that 
k n = k A K n , where K n is a vector of the reciprocal lat¬ 
tice.) However, since all functions (f> St i{k n , R) with dif¬ 
ferent index n have the same radial part, U s,1 (r) = u e A 
a s ,i u s A& s ,i ii s for i = 1, or U s ’ 2 (r) = ii e Aa s , 2 u s +b s , 2 u s 
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for i = 2, they are simply proportional to each other, 
4> s ,i{kni R) ~ ' t>s,i{k n 'i R)- Therefore, to avoid the linear 
dependence we should choose only one set of functions 
4>s,i(k n , R) corresponding to a single index n. The obvi¬ 
ous choice is to use the function 4> s ,i{k , R) with n = 0, 
ko = k and the reciprocal lattice vector Kq = 0. In that 
case the coefficient Co = n/V i l R 2 MT Y*m(k) [compare 
with Eq. (11611 ] can be further rationalized by omitting the 
multiplier Y* m (k) [or equivalently, including it in factors 
t a and tb, Eq. (Ifflap . (jlObjl ], Thus, we substitute Co with 

4tt 

Cq = -j=i R 2 mt . (19) 

(In principle, since the local function is not orthonormal, 
we can simply put Co = 1, but the form (1191) being similar 
to the constant coefficient for the standard LAPW basis 
function, simplifies some expressions for programming.) 

To study the transformational properties of supple¬ 
mentary basis functions 4> s ,i{k , R), Eq. (fT71) . we first 
rewrite them in the following form, 

R) = ti ,m (R - Ra), (20) 

where for each site a we have introduced two local func¬ 
tions (i = 1,2) of the Z, m-type, 

il>\’ m (R - Ra) = C 0 U s,i (r) Yr(f). (21) 

Notice that each local function ipl’ m (R — R a ) is strictly 
confined inside the MT-sphere a, because both U s,l {r) 
and Ri’ l (r) satisfy the boundary conditions (|18a[> . (|18bj) . 
We can then extend the function <f> s ,i{k , R) to the inter¬ 
stitial region ( R e I) by requiring 4> S: i(k , R) = 0. For 
the whole crystal we thus have 

Re MR ( 

R) = 0, Re I. 

This is a clear manifestation of the tight binding wave 
function. The multiplier e®*^“ in Eqs. (1221) and (12U1) en¬ 
sures that the supplementary wave functions 4> St i(k 1 R) 
obey the Bloch theorem. It is worth noting that usually 
the tight-binding description is spoiled by the presence of 
overlap between tails of wave function centered at neigh¬ 
boring sites. In the present method the tight-binding 
functions, Eqs. (1221) and ©, are free from this disadvan¬ 
tage because the local functions (and their first deriva¬ 
tives) go to zero at the sphere boundary and the overlap 
is absent. Thus, supplementary tight-binding functions 
can be considered as additional basis states orthogonal 
to the standard LAPW basis set. 

All matrix elements between the supplementary basis 
functions in the spherically symmetric potential are 
quoted explicitly in Appendix [B] and all matrix elements 
between (j) s ^ and standard LAPW basis functions (f> n are 
listed in Appendix [Cj For briefness we do not quote here 


the partial charges and electron density associated with 
supplementary basis states. They are tightly connected 
with the overlap matrix elements given by Eqs. (IB2al) . 
(IB4aD and (IB5al) of Appendix iBl and Eqs. (IC2al) . (IC3al) 
of Appendix [C] Concerning the full potential expressions 
for the extended basis set it is worth noting that after 
some algebraic transformations the equations can be ob¬ 
tained by selecting in standard FLAPW equations the 
contributions with the orbital indices l,m referring to 
the components of supplemented states and combining 
them together according to Eqs. (I15bl) . (I15cl) . (TT71) . 

The tight binding basis functions have a very impor¬ 
tant and practical property: they work even in the case 
when their expansion energy E s lies not far from the 
LAPW linear expansion energy E e . (We recall that the 
whole procedure is designed to treat the complicated case 
of semicore states when E s is supposed to be separated 
from E e by at least 10 eV.) The limiting case E s ss E e 
is considered in detail in Appendix [Dj and also discussed 
in calculations of Cd in Sec. nvn 


IV. PRACTICAL IMPLEMENTATION 
A. Face centered cubic structure of La 

We have applied the method developed in Sec. unto 
full potential electron band structure calculations of face 
centered cubic (fee) structure of lanthanum. Atomic lan¬ 
thanum has completely filled 5 p semicore electron shell 
lying at -22.12 eV which can slightly mix with valence 
states (5 d, 6s, 4/) at energies from -3 to -2 eV. For 
lanthanum here and below we use the Perdew-Burke- 
Ernzerhof (PBE) [l5| variant of the generalized gradient 
approximation (GGA) which gives rather accurate lattice 
constants for our systems. 

We have employed our original version of LAPW code 
with the potential of general form }lfil | . Integration in the 
irreducible part of the Brillouin zone has been performed 
over 240 special points. Angular expansions for the elec¬ 
tron density and wave function inside MT-sphere have 
been done up to L max = 8. The number of basis func¬ 
tions has been limited by the condition K max RMT = 9.0, 
resulting in 65 basis states. In addition to the standard 
LAPW basis functions we have considered 6 supplemen¬ 
tary tight binding basis functions with the p— angular 
dependence, which are located strictly inside the MT- 
sphere, Sec. m 

The 5 p semicore states have been treated as band 
states with the LAPW linear expansion energy E e lying 
0.5 eV above the p— band bottom energy (which is —9.5 
eV for the equilibrium lattice constant a = 5.315 A). The 
linear expansion energy for supplementary tight binding 
p— states has been fixed at 1 eV below the Fermi energy, 
E s (p) = 6.10 eV. Under these conditions two supplemen¬ 
tary radial functions shown in Fig. [T| are given by 

Ri(r) = u e (r) + a s pu s (r) + b s ^u s (r), (23a) 























6 



0 1,-2 3 r 4 

MT 


FIG. 1: Two supplemented tight binding radial functions of 
p —type for face centered cubic (fee) structure of lanthanum 
( a = 5.315 A). Radius is given in A, Tmt = 3.474 A stands 
for the muffin-tin radius. 

where a s q = 0.1940, b Si i = 0.4041, and 

i? 2 (0 = u e {r) + a s , 2 u s (r) + b s , 2 u s {r), (23b) 

where a S:2 = —1.9609, b St2 = 1.0705. Notice that the 
number of nodes for both radial functions is three (ex¬ 
cluding points with r = 0 and r = Rmt) which allows 
us to consider these functions as “compressed” 5 p basis 
states (i.e. with the principal quantum number n = 5) 
strictly confined within MT-sphere. The iteration proce¬ 
dure with supplemented tight binding functions has been 
stable converging to a self consistent solution without ad¬ 
ditional difficulties. Our PBE-GGA calculations result in 
equilibrium fee lattice constant a = 5.315 A which com¬ 
pares well with the experimental value, a exp = 5.304 A, 
Ref. [l3- 

To compare our treatment with the standard 
(LAPW+LO) method which uses only the first local 
function R\ (r) . Eq. (I23al) . we have performed a series 
of calculations, the results of which are summarized in 
Tables E El |ml In all cases the present method gives 
lower values of the total energy, Table Q] However, the 
energy difference which amounts to 0.145 eV for the poor 
basis set ( K max RMT = 7) becomes smaller for the inter¬ 
mediate basis sets and decreases to a small value (0.004 
eV) for the best basis set ( K max RMT = 9). Neverthe¬ 
less, even this difference is clearly noticeable in energy 


TABLE I: Total energy ( Etot , in eV) for various basis sets for 
fee calculations of La. a = 5.315 A, Rmt = 3.474 A, Eq = 
-231170 eV, A E = £ tot (FLAPW++) - £ tot (FLAPW+). 
FLAPW+-1- stands for the present scheme (FLAPW + 
2STBFs) with two radial functions and FLAPW+ for the 
FLAPW + LO method with a single radial function. 


Rmt ' Kmax 

FLAPW++ 

FLAPW+ 

A E 

7.0 

E 0 - 8.6878 

Eq - 8.5432 

-0.1446 

7.5 

E 0 - 9.2282 

Eq - 9.1730 

-0.0552 

8.0 

E 0 - 9.5362 

Eq - 9.5301 

-0.0061 

9.0 

E 0 - 9.5635 

Eq - 9.5599 

-0.0036 


TABLE II: Energy parameters (in eV) for various ba¬ 
sis sets for fee calculations of La (a = 5.315 A, Rmt = 
3.474 A). FLAPW++ stands for the present scheme (FLAPW 
+ 2STBFs) with two radial functions and FLAPW+ for the 
FLAPW + LO method with a single radial function. 



semicore 

valence 


5p—band 

(spd )—band 

Rmt Kmax 

Ebot Etop 

Ebot Ef 

7.0 

-11.1905 -9.9254 

2.7704 6.2047 

FLAPW++ 7.5 

-10.9033 -9.5956 

2.9965 6.3020 

8.0 

-10.0729 -8.7075 

3.7651 7.0381 

9.0 

-9.9972 -8.6297 

3.8310 7.1009 

7.0 

-10.9433 -9.6852 

2.9275 6.3934 

FLAPW+ 7.5 

-10.7326 -9.4289 

3.1170 6.4386 

8.0 

-9.9564 -8.5890 

3.8695 7.1478 

9.0 

-9.9240 -8.5550 

3.8958 7.1695 


band characteristics. In particular, band energy spec¬ 
trum demonstrates that the energy difference is of «0.07 
eV for the best basis set, Tables lIlT 1III1 It is worth noting 
that for all basis sets the present method results gives 
smaller band energies, Table [TT1 


TABLE III: Energy band spectrum (in eV) of fee La (a = 
5.315 A, Rmt = 3.474 A, Rmt ■ K m ax = 9) at the T-point 
of the Brillouin zone. FLAPW+-1- stands for the present 
scheme (FLAPW + 2STBFs) with two radial functions and 


FLAPW+ for the FLAPW + LO method with 
function. 

a single radial 

band 

deg. 

FLAPW++ 

FLAPW+ 

1 

(3) 

-8.6297 

-8.5550 

2 


3.8310 

3.8958 

E f 


7.1009 

7.1695 

3 


8.5444 

8.6258 

5 

(3) 

8.6128 

8.6803 




























7 




FIG. 2: Electronic band structure of bcc 7—La along high 
symmetry lines of the Brillouin zone (a = 4.243 A), (a) va¬ 
lence bands, (b) semicore (5 p) band. The horizontal dashed 
line indicates the Fermi level. 


B. Body centered cubic structure of La 

For the body centered cubic (bcc) phase of lan¬ 
thanum we have used a plane wave cut off parameter 
Kmax^MT = 9 (79 basis states plus 6 additional tight 
binding p— states), L max = 8 for the expansion of the 
electron density and wave functions inside MT-spheres, 
285 special points in the irreducible part of the Brillouin 
zone during the self-consistent procedure. 

As for fcc-La the LAPW linear expansion energy for 
the extended p— states of bcc-La has been has been cho¬ 
sen at 0.5 eV above the p— band bottom energy (i.e. 
E e (p) = —9.29 eV for the equilibrium lattice constant 
a = 4.243 A). The linear expansion energy for the sup¬ 
plemented tight binding p— states has been fixed at 1.0 eV 
below the Fermi energy ( E s (p ) = 6.32 eV). Two supple¬ 
mented radial functions, quoted in Eqs. (I23al) and (I23bl) . 
are defined by the coefficients a Sj i = 0.1502, b S: i = 0.4844 
for Ri(r), and a Si 2 = —1.8251, 6 Si 2 = 0.7718 for i? 2 (f)- 
The equilibrium lattice constant found for bcc-La, 
a = 4.243 A, is in good correspondence with the exper¬ 
imental value a exp = 4.25 A [18j. The calculated band 
structure of bcc lattice of lanthanum (PBE exchange and 
correlation) is plotted in Fig. [2] 

The results of the present and LAPW+LO approaches 
are compared in Tables EY1 El and eh As for fcc-La, 
the present scheme gives lower values of total energy of 


TABLE IV: Total energy ( Etot , in eV) for various basis sets 
for bcc calculations of La. a = 4.243 A, Rmt = 3.355 A, Eq = 
-231170 eV, A E = £ to t(FLAPW++) - £ tot (FLAPW+). 
FLAPW+-f- stands for the present scheme (FLAPW + 
2STBFs) with two radial functions and FLAPW+ for the 
FLAPW + LO method with a single radial function. 


Rmt ' Kmax 

FLAPW++ 

FLAPW+ 

A E 

7.0 

E 0 - 8.8166 

Eq - 8.7373 

-0.0793 

8.0 

E 0 - 9.3567 

Eq - 9.3428 

-0.0139 

9.0 

E 0 - 9.4574 

Eq - 9.4545 

-0.0029 


TABLE V: Energy parameters (in eV) for various basis 
sets for bcc calculations of La (a = 4.243 A, Rmt = 
3.355 A). FLAPW+-1- stands for the present scheme (FLAPW 
+ 2STBFs) with two radial functions and FLAPW+ for the 
FLAPW + LO method with a single radial function. 



semicore 

valence 


5p—band 

(spd)— band 

Rmt Kmax 

Ebot Etop 

Ebot Ef 

7.0 

-10.7891 -9.5738 

3.1271 6.5415 

FLAPW++ 8.0 

-10.2179 -8.9376 

3.6086 6.9193 

9.0 

-9.7907 -8.4968 

4.0213 7.3193 

7.0 

-10.5792 -9.3654 

3.2761 6.7062 

FLAPW+ 8.0 

-10.0659 -8.7832 

3.7403 7.0573 

9.0 

-9.7178 -8.4227 

4.0875 7.3877 


bcc-La for all basis sets. The total energy difference in¬ 
creases with worsening of the basis quality, Table IIVI For 
the best basis set ( Rmt • K max = 9) the present method 
gives energy spectrum shifted downwards by ~ 0.07 eV in 
comparison with FLAPW+LO band energy values, Ta¬ 
bles IVII El The difference between energy bands param¬ 
eters reaches 0.2 eV for poor basis set ( Rmt • K max = 7), 
Table El For the best basis set the total energy of bcc- 
La ( 7 —phase) is ~ 0.1 eV higher than the total energy of 
fcc-La (/3— phase), Tables HI and IIVI which is in agreement 
with the fact that at normal pressure 7 —La exists only 


TABLE VI: Energy band spectrum (in eV) of bcc La (a = 
4.243 A, Rmt = 3.355 A, Rmt ■ K max = 9) at the T-point 
of the Brillouin zone. FLAPW++ stands for the present 
scheme (FLAPW + 2STBFs) with two radial functions and 
FLAPW+ for the FLAPW + LO method with a single radial 


function. 

band 

deg. 

FLAPW++ 

FLAPW+ 

1 

(3) 

-8.4968 

-8.4227 

2 


4.0213 

4.0875 

E f 


7.3193 

7.3877 

3 

(3) 

8.8077 

8.8758 

5 

(3) 

8.9685 

9.0497 
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FIG. 3: (a) Electronic band structure of hexagonal close 

packed Cd along high symmetry lines of the Brillouin zone 
(a = 2.986 A, c= 5.632 A). 

at high temperatures (>1138 K) in narrow temperature 
range (53 K) [H]. 

C. Hexagonal close packed structure of Cd 

Hexagonal close packed (hep) cadmium is a special 
case because its Ad— states are not separated from the 
valence 5s—states by an energy gap, Fig. [3] In fact there 
is a small overlap between the top of d— band and the 
bottom of s—band. This implies that if the linear ex¬ 
pansion energy of extended d— states is properly chosen 
(for example, at 0.5 eV above the d— band bottom en¬ 
ergy, E(d) = —1.97 eV) the electronic band structure 
of cadmium can be carried out without supplemented 
tight-binding basis states of d— type. Thus, both calcu¬ 
lations i.e. with and without supplemented states, can 
be directly compared with each other. The second pe¬ 
culiarity is that because of the gapless energy spectrum 
the radial distribution of d— states does not change much 
throughout the valence band. If we chose the linear ex¬ 
pansion energy of supplemented d— states at 1.0 eV be¬ 
low the Fermi energy (E s (d) = 6.42 eV), then the en¬ 
ergy difference E s {d) — E e (d) ~ 8.5 eV is a relatively 
small value. In that case one can expect that the ex¬ 
tended d— states and the supplemented d— states are close 
to being linearly dependent, and we can test the sce¬ 
nario described in Appendix |Dj (Two supplementary ra¬ 
dial functions, Eqs. (123al) and (I23bl) . are specified by the 
coefficients a s ,i = —0.3448, 6 s p = 0.4643 for Ri(r) and 
a s>2 = -1.9343, b 8 , 2 = -0.2952 for R 2 (r).) 

The technical parameters of FLAPW calculations 
were the following: the plane wave cut off parameter 
KmaxRMT < 9 (149 basis states and 10 supplemented 
tight binding d— states), L max = 8 for the expansion 
of electron density and wave functions inside the MT- 
spheres, 216 special points in the irreducible part of the 
Brillouin zone during the self-consistent procedure, and 


TABLE VII: Total energy ( Etot , in eV) for various ba¬ 
sis sets for hep calculations of Cd. a = 2.986 A, c = 
5.632 A, Rmt = 2.74 A, E 0 = -304510 eV, A E = 
E tot (FLAPW++) - Etot (FLAPW+). FLAPW++ stands 
for the present scheme (FLAPW + 2STBFs) with two ra¬ 
dial functions and FLAPW+ for the FLAPW + LO method 
with a single radial function. 


Rmt Kmax 

FLAPW++ 

FLAPW+ 

A E 

FLAPW 

7.0 

Eo 

- 4.3780 

Eo 

- 3.8336 

-0.5444 

Eo 

- 2.4724 

7.5 

E 0 

- 6.6012 

Eo 

- 6.4998 

-0.1014 

Eo 

- 6.3540 

8.0 

Eo 

- 7.4471 

E 0 

- 7.4128 

-0.0343 

E 0 

- 7.3592 

8.5 

Eo 

- 7.9324 

Eo 

- 7.9033 

-0.0291 

Eo 

- 7.8428 

9.0 

Eo 

- 8.0093 

Eo 

- 7.9899 

-0.0194 

Eo 

- 7.9459 


TABLE VIII: Energy parameters (in eV) for various basis sets 
for hep calculations of Cd. FLAPW+-f- stands for the present 
scheme (FLAPW + 2STBFs) with two radial functions and 
FLAPW+ for the FLAPW + LO method with a single radial 
function. 


RmT ’ E-max 

Rbot 

Ef 


7.0 

-3.1660 

6.8094 


7.5 

-2.9380 

7.1192 

FLAPW++ 

8.0 

-2.7794 

7.2053 


8.5 

-2.5564 

7.3358 


9.0 

-2.4655 

7.4186 


7.0 

-2.9447 

6.9453 


7.5 

-2.8216 

7.2016 

FLAPW+ 

8.0 

-2.6999 

7.2648 


8.5 

-2.4763 

7.3969 


9.0 

-2.3949 

7.4741 


7.0 

-3.3958 

7.0200 


7.5 

-2.9187 

7.2652 

FLAPW 

8.0 

-2.7556 

7.2997 


8.5 

-2.5211 

7.4394 


9.0 

-2.4273 

7.5123 


the PBE m form of the exchange correlation potential. 

The equilibrium lattice constants are a = 2.986 A, c = 
5.632 A, which are in good correspondence with the ex¬ 
perimental values, a exp = 2.9794 A and c exp = 5.6186 A 
jl9l l. Band energies are shown in Fig. [3] 

Comparison between the present approach, the 
LAPW+LO and the standard LAPW treatment is pre¬ 
sented in Tables IVTIllVIHI and llXI We observe that for all 
basis sets the present approach gives lower total energy 
values, Table IVTII For the best basis set (RhirKmax = 9) 
the obtained total energy difference with the LAPW+LO 
value, 0.02 eV, is larger than for fcc-La or bcc-La. In 
the present approach band energies computed with the 
best basis set are lowered by 0.04-0.09 eV in comparison 
with LAPW+LO and LAPW values, Tables IVIIIl and UXl 
These energy shifts are also typical for the other points 
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TABLE IX: Energy band spectrum (in eV) of hep struc¬ 
ture of Cd at the T—point of the Brillouin zone (o = 2.986 A, 
c = 5.632 A) with PBE exchange 0 ( Ef is the Fermi energy, 
deg. is the energy degeneracy). FLAPW+-(- stands for the 
present scheme (FLAPW + 2STBFs) with two radial func¬ 
tions and FLAPW+ for the FLAPW + LO method with a 
single radial function. 



deg. 

FLAPW++ 

FLAPW+ 

FLAPW 

1 


-2.3356 

-2.2662 

-2.2954 

2 


-2.0708 

-2.0009 

-2.0305 

3 

(2) 

-1.6919 

-1.6194 

-1.6604 

4 

(2) 

-1.6399 

-1.5668 

-1.6078 

5 


-1.4047 

-1.3485 

-1.3215 

6 

(2) 

-1.1638 

-1.0888 

-1.1334 

7 

(2) 

-0.9125 

-0.8377 

-0.8866 

8 


4.3215 

4.3750 

4.4176 

9 


4.4792 

4.5374 

4.5479 

Ef 


7.4186 

7.4741 

7.5123 

10 


17.6841 

17.7375 

17.7711 


of the Brillouin zone. 


V. CONCLUSIONS 

We have presented a new method for the improvement 
of the LAPW description of the electronic band structure 
by using two linearization energies for the same (l,m) 
partial component. Starting with two LAPW radial func¬ 
tions, Eqs. (I5bl) and (l5cl) . having the same angular depen¬ 
dence Yi ctm (r) but different linearization energies (E^ 
and Ej' 2> ) inside MT-spheres, we have demonstrated that 
their augmentation to the basis plane wave can be per¬ 
formed by constructing additional basis functions <j> 8t i 
(i = 1,2) in the form of Eq. (fT71) [two functions, Eqs. 
(I15bl) and (I15cl) . for each l, m-component]. The sup¬ 
plementary basis functions have zero values and slopes 
on the sphere surface, Eqs. (118al) . (118b[) . and are linear 
independent of the usual LAPW basis states. The con¬ 
structed basis functions are of the tight-binding type, Eq. 
(l22l) . and obey Bloch’s law. 

In contrast to the LAPW+LO method with only one 
supplemented function, Eq. (I15bl) . in our treatment for 
each l, m-component there are two supplemented func¬ 
tions [Eq. (I15bl) and (I15c[) ]. The second basis function 
(absent in LAPW+LO) closely examined in this work, 
owes its appearance to the ui function in the canoni¬ 
cal LAPW method. Thus, the basis sets of LAPW and 
LAPW+LO methods can be extended further by adding 
supplemented functions of the tight-binding type, Eq. 
(I15cl) . 

In Sec. lIVl the present method with extended basis set 
has been applied to the study of the face centered and 
body centered phases of lanthanum (/?—La and 7 —La) 


with the 5p— semicore shell separated by a gap of for¬ 
bidden states from the valence states and to the hexago¬ 
nal close packed structure of cadmium, where the semi¬ 
core 4d—states overlap with the valence 5s—states. In 
all cases we have observed a systematic improvement 
in the values of total energy in comparison with the 
standard LAPW+LO treatment, Tables ffluaEm The 
difference with LAPW+LO total energy is only 0.003- 
0.004 eV for La and 0.019 eV for Cd for the best basis 
set ( Rmt • K max = 9) but significantly increases in go¬ 
ing to intermediate ( Rmt ■ K max = 8.5 or 8 ) and poor 
{Rmt ■ K max = 7.5 or 7) basis sets. 
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Appendix A 


Solution to the system of linear equations, Eq. 
is given by 

pii, 

0+1 = Ueil's - u e u s ), 

(Ala) 

b 8 ,1 = u' e u s - u e u' s ), 

(Alb) 

and the solution to the svstem (Il2bl) is 


a s ,2 = {u e u' s - u e u s ), 

(A2a) 

b s ,2 = ^(«' e u s - u e u' s ). 

(A2b) 

Here 


A • / •/ 1 

A = u s u s u s u s « 2 . 

K MT 

(A3) 

Appendix B 


The matrix elements for the overlap and Hamiltonian 
operator between supplementary states, 

= oii, 

(<)>.,i\H\</>s tj ) = HH, 

(Bla) 

(Bib) 


are partitioned in three different blocks, when i = j = 1 
(block I), i = j = 2 (block II), and i = 1, j = 2 or i = 2, 
j = 1 (block III). 

For the first block (/) we have 

Oll(k) = O s , s (1 + c£i + 0, 


(B2a) 
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where N a is the number of equivalent spheres a, while 

O s>s = N a , (B2b) 

C s s — T b s ,ib s ^i J\f (u s , u s ), (B2c) 

Ces — ^s,lA/" (Ue, u s ) T b s ^l J\f (u e , U s ), (B2d) 

and M stands for the integral over the product of two 
functions, 

r^MT 

N(u\,u 2 ) = / ui(r)u 2 (r)r 2 dr. (B3) 

Jo 

Notice that N(u e ,u e ) = A f(u s ,u s ) = 1 and A f(u e ,u e ) = 
J\f(u s ,u s ) = 0. It is also assumed here that the constant 
coefficient of the supplementary function is taken in the 
form of Eq. (fH)l) . If another form is used, the factor 0 StS , 
Eq. (IB2bl) . [and O e ^(k n ), Eq. (IC'2cl) ] should be changed 
accordingly. 

For the block II we obtain 



Of s = 0 S , S (A f(u e , lie) + C 22 s + Cel), 

(B4a) 

where 


C 2 s 2 s = a s , 2 a s , 2 + bs, 2 b s , 2 J\f(u s ,Us), 

(B4b) 


Ce 2 s = a s , 2 J\f(u e ,u s ) + b St 2 J\f(u e , u s )- 

(B4c) 

Finally, for the block III we get 



Oil = Os,s (Cll + Cel + Cel), 

(B5a) 

where 


^ss — ^S,1^S,2 ^S,1^S,2 -hf (ils ? ^s) ? 

(B5b) 


and C e l i s given by Eq. (IB2dl) . while C e 2 s by Eq. (IB4cl) . 

For the matrix elements of the Hamiltonian H l s i(k), 
Eq. (IBlbl) . we also obtain three blocks. For the first 
block (I) we have 

HU = 0 StS ( E e + E s CH + a Si i& Sj i 
+ (Ee + E s ) C e l + b St iN(u e , u s ) ) . (B6) 

Here E e and E s are energies at which the radial wave 
functions u e (r) and u s (r ) are evaluated in the MT -sphere 

a. 

For the second block (II) we get 
H 22 s = O a , a (. E e M(u e , ii e ) + E s Cf s + a s , 2 bs ,2 

+ (E e + E s ) Ce 2 s + 7 es) ) (B7a) 

where 

7es = «s ,2 A/"(w e , U s ) +b St 2 (N(u e ,U s ) + J\f(u e ,U s )) . 

(B7b) 

For the third block (III) we have 

Hl 2 s =O s ,s (E s (Cl 2 s + C e l) + a St 2 b s ,i 

+ E e C e l + A u s )l . (B8) 


Appendix C 

In this section we quote the expressions for matrix ele¬ 
ments for the overlap and Hamiltonian operator between 
supplementary and extended states, 

(<t> s ,i\0\(j) n ) = Oi’ n e(k), (Cla) 

(<t>s,i\H\K)=Hl’ n e (k). (Clb) 

The extended states here are the usual LAPW basis func¬ 
tions, Eq. HD, which are characterized by the wave vec¬ 
tor k n = k + K n . The supplementary functions have two 
components i = 1 , 2 , for each Z, m-angular dependence, 
Eq. (I15bl) . (I15cl) . 1171) . 

For the matrix of overlap we get 

Ol:(k)=O e , s (k n )Sl ”, (C2a) 

where 

Si™ = a e + a s ,ia e J\f(u e ,Us) + a Sj ib e M(u e , u s ) 

+b s ,ia e M(u ei iis) + b e bs,iJ\f(u e ,u s ), (C2b) 

and the structure factor is 

OeAkn) = ^(R a M T ) A Y l * m (k n ), 

x ^2 exp[i(k n - k) f v>a \. (C2c) 

V 

Here, v va stands for the coordinates of all v centers 
of MT-spheres of the type a in the primitive unit cell. 
The factors J\f(ui,u 2 ) in (IC2bl) are integrals between two 
functions given by Eq. dEap . 

The quantities a e and b e in Eq. (IC2bl) and below are 
standard LAPW expansion coefficients for the compo¬ 
nent with l , m, defined by Eq. (l9al) and (l9bP . i.e. a e = a° 
and b e = b° e . Since they depend on k n , l, and a, we can 
write a e = a^ik), b e = b^fk). Explicit expressions for 
them can be found in Ref. @1 ■ 

The matrix element for the second case (i = 2) reads 
as 

0 2 s n e (k)=0e,s(k n )S 2n e, (C3a) 

where 

S s e b e N(u e , Ue) T U Si2 U e A/* (u e , U s ) T Us^ 2 b e A/* (u e , U s ) 

+bs, 2 a e JV(ue,Us) + b e bs, 2 N(u e ,u s )- (C3b) 

Below we quote the matrix elements for the Hamilto¬ 
nian, 

Hle(k) = O e ,s(kn) (E e S]” + b e + a s ,ib e J\f(u e ,Us) 

+b s ,ib e J\f(u e ,Us)), (C4) 

H 2 s n e (k) = O e ,s(k n ) (Ee S 2 : + a s , 2 beM(Ue,Us) 

+bs t2 b e JV(u e ,Us))- (C5) 

Here O e , s is given by Eq. (IC2cl) . while S'*” by Eq. (IC2bl) 
for i = 1, and by Eq. (lC3bl) for i = 2. 
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Appendix D 

Here we demonstrate that the two supplementary basis 
functions <t> s ,i (i = 1,2), Eq. (fT71) . work even in the case 
when the expansion energies E s and E e lie not far from 
each other. 

Consider E e = E s + e, where e/E s <C 1. Making use 
of the following expansions 

u e (r) = u s (r) + u s (r) e 

+ i?i s (r) e 1 + ^ u s (r) e 3 + 0(e 4 5 6 7 8 9 10 11 ), (Dla) 

u e (r) = u s (r) + u s (r) e + ^u s (r) e 2 + 0(e 3 ), 

(Dlb) 

for u e (r) and u e (r) and substituting them in Eqs. (115bl) . 
(115cll . we arrive at 

R i’L(. r ) = \ £2 °o e ikRa [« s + u a + 0(e 2 ) 

+a' sl u s + b' sl u s \, (D2a) 

R i’m(. r ) = £ Co [ii s + — e ii s + 0(e 2 ) 

+a' s2 u s + b' s2 u s \. (D2b) 

Here 

a's,i = + !)> b 's,i = ( D3a ) 

a 's 2 = — a s,2, b' s 2 = - (6 s ,2 + 1). (D3b) 

E ’ £ 

The prefactors e 2 /2 and e in Eqs. (ID2al) . (ID2bl) are not 
very important, because as a consequence of solving sec¬ 


ular equations these basis functions will be effectively or- 
thonormalized. Functions ^’^(r) and R{’ m (r) have two 
important features. First, by a linear transformation the 
two functions, Eqs. (ID2al) . (ID2bl) . can be transformed to 
two functions with the following radial dependencies 

U\ (r) = u s + a" 1 u s + b" 1 u s , (D4a) 

U 2 (r) = «, + a" tl u s + 6 "! u s , (D4b) 

where the coefficients a" x etc. can be expressed through 
a' s 1 etc. Since the functions u s (r) and w s (r) are linear in¬ 
dependent, the same property applies to U\{r) and Uz(r) 
and, consequently, to l ( r ) and 7?, / s m (r), although the 
property deteriorates as e —> 0. Second, we still can im¬ 
pose the boundary conditions (118al) and (118bll . 

However, since the initial functions i?f’ TO ( r ) and 

R i m ( r ), have almost identical radial dependence (ne¬ 
glecting terms with e u s and the others of high order of 
e in Eqs. (ID2al) . (lD2bl) make them completely identical), 
the normalization procedure leads to the appearance of 
an effective basis state </>', which is orthogonal to other 
basis states and expressed through the linear combina¬ 
tion of initial states, <j)' = C\ ip St i + C2 ips,2 with large 
coefficients C\ and C2 (i.e. \C\\ 1, \Cy,\ 1). In that 

case, the partial charges of supplementary basis states 
can also be very large. Nevertheless, some of these par¬ 
tial charges are of opposite sign and effectively cancel 
each other in the final answer. We have observed the 
effect in the calculation of the hexagonal close packed 
lattice of cadmium reported in Sec. II V Cl 
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